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REDUCED COMPLEXITY MLSE EQUALIZER 
FOR M ARY MODULATED SIGNALS 

BACKGROUND 

Digital wireless communication systems are used to convey a variety of 
information between multiple locations. With digital communications, information is 
translated into a digital or binary form, referred to as bits, for communication 
purposes. The transmitter maps this bit stream into a modulated symbol stream, which 
is detected at the digital receiver and mapped back into bits and information. 

In digital communications, the radio environment presents many difficulties that 
impede successful communications. One difficulty is that the signal level can fade 
because the signal may travel in multiple paths. As a result, signal images may arrive 
at the receiver antenna out of phase. This type of fading is commonly referred to as 
Rayleigh fading or fast fading. When the signal fades, the signal-to-noise ratio 
becomes lower, causing degradation in the quality of the communication link. 

A second problem occurs when the multiple signal paths are much different in 
length. In this case, time dispersion occurs, in which multiple fading signal images 
arrive at the receiver antenna at different times, thus giving rise to signal echoes or 
rays. This causes intersymbol interference (ISI), where the echoes of one symbol 
interfere with subsequent symbols. 

At the receiver, coherent demodulation is desirable, since it provides the best 
performance. This requires knowledge of the multipath channel. In many wireless 
applications, this channel is time-varying, due to transmitter motion, receiver motion, 
and/or scatterer motion. Thus, there is a need to track a time varying multipath 
channel. 

To provide coherent demodulation of multipath signals, a maximum-likelihood- 
sequence-estimation (MLSE) equalizer may be employed. Such an equalizer considers 
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various hypotheses for the transmitted symbol sequence, and with a model of the 
dispersive channel, determines which hypothesis best fits the received data. This is 
efficiently realized using the Viterbi algorithm. This equalization technique is well 
known to those skilled in the art, and can be found in J.C. Proakis, Digital 
5 Communications, 1989. 

The conventional MLSE equalizer can be explained by a simple example. 
Suppose the transmitter transmits a symbol stream s(ri), which takes on values " +B" or 
"-2T corresponding to bit values 0 or 1, respectively. This stream is modulated using 
binary-shift keying (BPSK). At the receiver, the received signal is filtered, amplified, 
10 and mixed down using the I and Q carriers, then sampled once every symbol period 
(7), giving a received signal stream r(ri). In this example, the intervening channel 
consists of two rays, a main ray and an echo, where the echo arrives T seconds later 
and 7* is the symbol period. Then the received signal can be modeled as 

r(ri) = CosQi) + c x s(n - 1) +r|(n), (1) 

15 where c 0 and c x are complex channel tap values and r|(n) is additive noise or 
interference. 

In the MLSE equalizer, at iteration n, there would be two different previous 
"states," 0 and 1, corresponding to the two possible values for the previous symbols: 
1. s(n - 1) = B - Previous State = 0 
20 2. s(n - 1) = -B « Previous State = 1. 

Associated with each previous state there would be an accumulated metric, 
accumulated from previous iterations, giving rise to accumulated metrics A 0 (n - 1) and 
A x {n~ 1) for previous state 0 and previous state 1, respectively. 

There would also be two current states corresponding to two possible values for 
25 s{n). Each possible pairing of a previous state with a current state corresponds to a 
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hypothetical sequence {s h (n - 1), s h (ri)}. For each of these hypotheses, there will be a 
corresponding hypothesized received signal value at time n: 

r h {n) = c& h (n) + c x s h (n - 1). (2) 

Furthermore, for each of these hypotheses, there will be a corresponding "branch" 
metric given by 

M h {n) = \r(n)~r h (n)\\ (3) 

In the example, there are four possible hypothetical sequences which can be 
denoted by h e {00, 01, 10, 11}, as illustrated in Table L 



Previous State 


Current State 


h 


0(B) 


0(B) 


00 


1(-B) 


0(B) 


10 


0(B) 


K-B) 


01 


K-B) 


K-B) 


11 



Table 1: Naming Convention for States 



The candidate metrics for each possible current state would be the sum of the 
corresponding branch metrics and the previously accumulated metric associated with 
s h (n - 1). For each current state, there are two possible previous states. For each 
current state, the previous state which gives the smallest candidate metric is selected as 
the predecessor state, and the candidate metric becomes the accumulated metric for that 
current state. 
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Thus, for current state 0 (i.e., s(n)=B), there will be two hypothetical 
sequences {s h (n - 1) = B, s h (n) = B) or {s h (n - 1) = -B, s h (n) = B}, denoted by 
hypotheses h = 00 and h = 10. This gives rise to two candidate metrics for current 
state 0: 



5 Coo(tt) = A 0 (n - 1) + M m (n) (4) 

C 10 (n) = A x {n - 1) + M 10 (n). (5) 

The smaller of these two candidate metrics gives the accumulated metric for 
current state 0. The corresponding previous state becomes the predecessor state to 
state 0 at time n: 
10 1. If C 10 (n) < C m {n),A^ri) - C 10 (n) 

2. If Coo(n) < C 10 (n), A 0 (n) = C^n) 

A similar procedure is applied to current state 1 (i.e., s(n) = -B), in which case 
the two hypotheses are h = 01 and h = 11. 

For each n, four hypothesized received signal values must be computed 
15 according to equation (2). Computing each hypothesized value r h (n) involves two 

multiplication operations. Each of these multiplication operations involves the product 
of one of the channel taps and a hypothetical symbol value. 

If B is set to + 1, the multiplications become simply sign changes. However, 
if amplitude modulation is used, this is not possible. Also, suppose symbol values can 
20 be normalized to j (i.e., 8PSK modulation). In this case, the necessary 

multiplications can not be^ implemented as simple sign changes. 

One way to increase the throughput over a digital channel is to allow each 
transmitted symbol to take on more than two values. In general, each symbol might be 
chosen from a set of M possible values {B 1? B 2 , . . . , B M } thereby increasing the number 
25 of bits per symbol. The complexity of the conventional implementation of the MLSE 
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equalizer increases exponentially with the number of allowed symbols (i.e., with M). 
Specifically, with "L" channel taps, the number of multiplications needed for the 
conventional implementation of the equalizer is proportional (proportionality constant 
greater than one) to M w , 

5 Even with moderate values for M and L (e.g., 8 and 5 respectively), it is 

impractical to implement the equalizer in grand purpose DSPs with low power 
consumption. The trend in power sensitive applications is to move away from DSPs 
toward ASICs. ASICs can be optimized for performing one task with low power 
consumption and small size (i.e., small gate count). Implementing multipliers in an 

10 ASIC is much more expensive (in terms of power consumption and gate count) than 
implementing circuits for addition or subtraction. Therefore, it is highly desirable to 
avoid many of the multiplications associated with the conventional implementation of 
MLSE equalizers. 

Because of the trend toward increasing the number of required multiplication 
15 operations and due to the fact that technology is moving away from DSPs, where 

multiplication operations are easy, toward ASICs, where multiplications are relatively 
expensive, it is desirable to reduce the number of multiplications actually performed by 
the MLSE. 



SUMMARY 

20 This disclosure concerns demodulation of radio signals modulated with M-ary 

modulation in the presence of intersymbol interference distortion. The invention 
presents a method for reducing the number of multiplications needed to implement an 
MLSE equalizer for signals modulated with M-ary modulation. 

In exemplary embodiments of the present invention, the number of 

25 multiplications is reduced by pre-computing certain values needed for the determination 
of the branch metric and storing these pre-computed values in a product table. When a 
branch metric computation is to be made, whether it is an Euclidean branch metric 
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computation or an Ungerboeck branch metric computation, certain multiplication 
operations are replaced by simple table look-up operations. As a result, the power 
consumption and size of the equalizer are reduced. 

Any receiver that demodulates signals that are modulated with M-ary 
5 modulation can be implemented using the instant invention. The resulting demodulator 
will have a lower complexity than existing demodulators. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The above objects and features of the present invention will be more apparent 
from the following description of the preferred embodiments with reference to the 
10 accompanying drawings, wherein: 

Figure 1 illustrates a radio communication system within which the techniques 
according to the present invention can be implemented; 

Figure 2 illustrates a functional block diagram of the baseband processor of 
Figure 1; 

15 Figure 3 illustrates a functional block diagram of the transmission function of 

Figure 1; 

Figure 4 illustrates a system for computing an Euclidean branch metric within 
which the techniques according to the instant invention may be implemented; 

Figure 5 illustrates a functional block diagram of a filter used to compute an 
20 Euclidean metric according to the prior art; 

Figure 6 illustrates a product table for the Euclidean metric according to one 
embodiment of the present invention; 

Figure 7 illustrates a functional block diagram of the filter for computing the 
Euclidean metric according to one embodiment of the present invention; 
25 Figure 8 illustrates a system for computing an Ungerboeck branch metric within 

which the techniques according to the instant invention may be implemented; 
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Figure 9 illustrates a functional block diagram of a conventional filter used to 
compute the Ungerboeck metric according to the prior art; 

Figure 10 illustrates a product table for the Ungerboeck metric according to one 
embodiment of the present invention; and 
5 Figure 1 1 illustrates a functional block diagram of the filter for computing the 

Ungerboeck metric according to one embodiment of the present invention. 

DETAILED DESCRIPTION 

ijj Figure 1 illustrates a radio communication system within which the present 

I*t invention may be implemented. In Figure 1, a radio transmitter and receiver for a 

10 radio communication system are provided. The radio communication system may 
IJj operate using FDMA, TDMA, or CDMA, or some combination thereof. A transmitter 

J has a digital symbol generator 102 which receives an information carrying signal 101 

H and generates a corresponding digital symbol sequence, S. The symbols S are 

Rf subjected to digital-to-analog (D/A) conversion, modulation, pulse shape filtering and 

^ 15 amplification, and are transmitted as analog signal Fby digital transmitter 103 
according to known techniques. 

Signal 7 travels through the radio channel and is intercepted by the antenna 104 
at the receiver. Thermal noise n is also intercepted by the antenna 104. 

Radio unit 105 amplifies, down-converts, and filters the received signal 
20 according to known methods to produce an analog output. This analog output is 
coupled to an A/D converter 106 which converts the analog signal into a received 
signal sample stream r(kT s ), where T s is the sample period, and it is an integer counter. 
The sampling period T s may be less than the symbol period T. The received signal 
sample streams are collected in processor 107, which processes this stream to produce 
25 an estimate of the transmitted digital symbol stream §. In later descriptions, 

transmission function 109 is used to refer to the signal path through digital transmitter 
103, the radio transmission channel 105, and A/D 106 collectively. 
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The transmission function 109 produces the received signal sample stream r(kT s ) 
which is sent to processing unit 107 where it is processed in accordance with the 
present invention. 

A functional block diagram of the baseband processing unit 107 is illustrated in 
5 Figure 2. The received signal sample stream r(kT s ) is coupled to a signal pre- 
processor, or sync, 206 where the received signal sample stream is correlated with a 
known timing/synchronization sequence according to known techniques. For the case 
of symbol-spaced demodulation, if the sample period T s is less than the symbol period 
y3 T, the signal pre-processor 206 performs a decimation of the received signal sample 

|H 10 stream r(kT s ) to produce one sample per symbol, designated as r(n). For fractionally - 

spaced demodulation, more than one sample per symbol is generated. 
«1 Estimating circuit 202 produces channel tap estimates c(z) which are used to 

4:1 * 

s model the radio transmission channel according to known techniques. This might 

2 involve an initial channel estimation period followed by a tracking period. The channel 

15 tap estimates c(t) are coupled to the input of the branch metric processor 203. The 
branch metric processor 203 is coupled to a sequence estimation processor 204 which 
w provides an estimate of the digital symbol stream S. 

The transmission function 109 is illustrated in more detail in Figure 3, where 
for simplicity, the number of received antennas is restricted to one. One skilled in the 
20 art will appreciate that the present invention may also be used for the case where there 
are two or more antennas. In Figure 3, the symbol sequence S is input to the digital 
transmitter 103 which transmits analog signal Y. The analog signal Y propagates 
through a radio transmission channel to the radio unit 105. Radio channel 301 may 
introduce fading and time dispersion. Omnipresent thermal noise n is also received. 
25 Radio unit 105 amplifies, down-converts, and filters the received signal according to 
known techniques to produce an analog signal. This analog signal is coupled to an 
A/D 106 which converts this analog signal into the received signal samples r(kT s ). 
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In an MLSE equalizer, all possible transmitted symbol sequences S are 
effectively considered. In one implementation, hypothesized symbol values s h (n) are 
filtered by channel tap estimates c(t) to produce hypothesized received samples r h (n). 
The difference between the hypothesized r h {n) and the actual r(n) received signal 
5 sample stream, referred to as the hypothesis error e(n), gives an indication of how good 
a particular hypothesis is. The squared magnitude of the hypothesis error is used as a 
metric to evaluate a particular hypothesis. The metric is accumulated for different 
hypotheses for use in determining which hypotheses are better using the sequence 
estimation algorithm. This process may be efficiently realized using the Viterbi 

10 algorithm. A description of the Viterbi algorithm can be found in G. Forney, Jr., "The 
Virterbi Algorithm," Proceedings of the IEEE, vol. 61, no. 3, March 1973, pp. 267- 
278. As will be appreciated by one skilled in the art, other sequence estimation 
algorithms may also be used. 

In an MLSE equalizer, there are states associated with different transmitted 

15 sequence hypotheses s h (ri). At a given iteration, there are previous states, each of 

which is associated with an accumulated metric. Each pairing of a previous state with 
a current state results in a branch metric M h (ri). The candidate metric for a current 
state is then the sum of this branch metric M h (n) and the previously accumulated 
metric. For each current state, the previous state which gives the smallest candidate 

20 metric is selected as the predecessor state, and the smallest candidate metric becomes 
the accumulated metric for the current state. The branch metric can be expressed as: 

M h {n)= \r(n) - r h (n)\ 2 (6) 

where 

jfc=0 
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The channel tap estimates are designated by c(t) where % is the delay (i.e. r = 0 is the 
main ray, x = 1 is the first echo, etc). N t is the number of channel taps estimated. For 
each n, the hypothesized received signal r h (n) must be computed according to equation 
(7). In the most general case, computing each hypothesized received value r h (n) 
involves N t complex multiplications (i.e. 4N t real multiplications). Each of these 
complex- valued multiplications involves the product of one of the estimated channel 
taps and a hypothesized transmitted signal. 

Figure 4 illustrates a system for computing an Euclidean branch metric M h (n) 
within which the techniques according to the instant invention may be implemented. In 
Figure 4, the hypothesized sequence of symbols s h (n) which is generated by symbol 
sequence generator 410 is coupled into filter 400 to produce the hypothesized received 
samples r h (n). The difference between the hypothesized received samples r h (n) and the 
actual received signal sample stream r(n) is the hypothesized error e(ri). The squared 
magnitude of the hypothesized error is performed by unit 403 to produce the branch 
metric M h (n). 

A functional block diagram of a conventional filter for computing an Euclidean 
branch metric is illustrated in Figure 5. In Figure 5, N t complex multiplications are 
performed to compute each r h (ri). If the transmitted symbols are restricted to be in a 
set of M possible values, {B l , 5 2 , . . . , B M } , the present invention provides a method for 
implementing the filter 400 that avoids performing any multiplications. 

With an Euclidean metric, assume that each hypothesized transmitted symbol 
s h (ri) is in the set {B x , B 2 , . . . , B M } . The first term in the sum in equation (7) is in the set 
{ B c(0)} M • A11 the members of this set can be pre-computed and stored in the first 
column of an M X N t table (see Figure 6). 

The j-th column of this table, corresponding to estimated channel tap c(j - 1), 
stores all the possible values of (s h (.)c(j - 1)), i.e. {B u c(j - 1), B 2 c(j - 1),..., B^fj - 1)}. 
Each hypothesized received value r h {n) can then be computed by simply adding the 
appropriate entries from this product table (see Figure 7). 
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As a concrete example, consider a two-tap channel (c 0 , c x ) with 8PSK 



modulation i.e., s h (n)e { /-^ 



2n 



«V . In 



this case, the first column of the table illustrated 



in Figure 6 will have 8 entries, corresponding to channel tap c 0 , as: 



.2jc rt 2n 2ti 2te 

j — 0 j j 2 / 7 

' 8 8 8 '8 

0' 0' 0' 5 0 



(8) 



Similarly, the second column of the table in Figure 6 will have 8 entries, corresponding 
to channel tap c u as: 



.27U 271 2Tt 2TC 

y 0 J J 2 / 7 

8 8 8 J 8 



(9) 



Symmetry can be used to reduce storage requirements. First, negative symmetry can 
be used to halve the items stored. Basically, 



2TC, 



27t„ 

e c 



(10) 



Therefore, only values for H — 0, 1, 2, 3 need to be stored. Furthermore, 



2TT 271 

/ 2 j 0 

8 . 8 

e c n = j e c 



(11) 



and 



e c = / e c 



(12) 



Thus, by switching real and imaginary parts and negating the new real part, one only 
needs to store values for « = 0 and 1. For H = 0, no multiplication is necessary since 
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2K 



e 




(13) 



Thus, one can simply store c 0 and e jn/4 c 0 and derive the other values using logic that 
negates and possibly switches real and imaginary parts. This savings results from the 
quadrantile symmetry in the 8PSK signal constellation. 

The Ungerboeck metric is obtained from the Euclidean metric M h (ri) of equation 
(3) in two steps. The first step is to expand M h (n) to get 

M h (n) = A(n) + B(n) + C(n) + D(n), (14) 

where 

A(n) = \r(n)\ 2 (15) 
B(n) = -2Re{r(n)cls* h (n)} -2Re{r(n)c *s * (n~l)} (16) 

C(n) = \c 0 \ 2 \s h (n)\ 2 + \ Cl \ 2 \ Sh (n-l)\ 2 (17) 
D(n)=2Re{c 0 cl s k (n)sl(n-l)}. (18) 



The Ungerboeck method drops the term A(n), which is common to all hypotheses. The 
second step combines terms proportional to s * h (n) from different iterations. At 
iteration (n + 1), the terms become: 

B(n + l) = -2Re{r(n + l)c* 0 s* (n +1)} -2Re{r(n+\)c *s * (n)} (19) 

C(n+1) = |c 0 | 2 |.y A («+l)| 2 + Ic^l^n)! 2 (20) 

D(n+ 1) = 2/?e{c 0 c; s h (n+ l)s* h («)}. (21) 
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Thus, there are terms proportional to s* h (n) in both iterations. These terms can 
be combined by defining a new metric, M h (n) , as 



M h (n) = Re {si (n) [-2z(n) + qo s h (n) + 4 A(n-l)]} (22) 

where 

Z(n) = c* Q r(n) + c\r{n + 1) (23) 

q 0 = |c 0 | 2 + |c,| 2 (24) 

o * (25) 
= 2c o c i ■ 

Here, q 0 and q! are referred to as s-parameters. 

Thus with an Ungerboeck metric, the branch metric M' h (n) is defined as: 

M'An) = Re{ S ;(n)[-2z(n)+t(n)]}, (26) 



where 



t h (n) = Y,^(k)s h {n-k). (27) 



Figure 8 illustrates a system for computing an Ungerboeck branch metric within 
which the techniques according to the instant invention may be implemented. In Figure 
8, the hypothesized symbols s h {n) are coupled into the filter 600, with impulse response 
q(r), to produce the t A (n). q(t)'s are determined from the estimated channel taps c(z)'s; 
hence, g(r)'s do not depend on the hypothesized symbols s h (n). The current state of the 
art for implementing a filter of Figure 8 requires N t complex multiplications, as 
illustrated in Figure 9. 
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Assume that each hypothesized transmitted symbol s k (n) is in the set {B x , B 2 ,..., 
B M }. The first term in the sum in equation (27) is in the set {B.qiO)}^ All the 
members of this set can be pre-computed and stored in the first column of an M X N t 
table (see Figure 10). 

5 The y-th column of this table, corresponding to q(j - 1), stores all the possible 

values of s h (.)q(j - 1), i.e., {B x q(j - 1), B 2 q(j - 1),..., B^qij - 1)}. Each hypothesized 
t h {ri) can then be computed by simply adding the appropriate entries from this product 
table (see Figure 11). 

As described earlier, quadrantile symmetry can be used to reduce the number of 

10 items stored in each column to two. 

The partial Ungerboeck metric also has terms like q t set forth in the previous 
section; hence, the present invention can be used to reduce the number of 
multiplications needed to implement a demodulator that uses the partial Ungerboeck 
metric. The partial Ungerboeck metric is disclosed in U.S. Patent 5,499,272 to 

15 Bottomley. 

For a multi-channel MLSE with Euclidean metric, the transmitted symbol 
stream is received via several channels. With two channels, there would be two 
received signals rjn) and rjn), corresponding to channels a and b respectively: 

rjn) = c a0 s(n) + c al s(n - I) + t]Jn) (28) 
20 r b (n) = c b0 s(n) -h c bI s(n - 1) + rj b (n). (29) 

Corresponding to each received signal and each hypothesized symbol sequence, 
there would be a branch error signal: 

ejn) = rjn) - r ah (n) (30) 
«wW = r b (n) -r bh (n), (31) 
25 where r ah (n) and r bh (n) are the hypothesized received signals for channel a and channel 
b (both corresponding to the same hypothesized sequence {s h {n - 1), s h {n - 1)}: 

rjn) - c<fi h (n) + C aI s h (n - 1) (32) 
r bh (n) = Cufi h (n) + C bl s h (n - 1). (33) 
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Note that one branch metric for channel a is e dh (n)e* ai (n) i and one branch metric for 
channel b is e bh (n)e* bh (n). 

Different branch metrics for the joint equalization of the outputs of both 
channels can be obtained by combining e ah (n) and e bh (n) in different ways. Several 
interesting branch metrics for the joint equalization problem can be expressed in the 
following form: 

M h (n) = E h (n) H QE h (n), (34) 

where 



EM = 



e ah( n ) 
e bh< n ) 



(35) 



and Q is a two by two weighting matrix. 
10 Any branch metric of the form given by equation (34) can be computed 

efficiently using two pre-computed tables. Each hypothesized received signal r ah (n) or 
r bh (n) can be computed without performing a multiplication operation by using two pre- 
computed tables. 

It would be appreciated by one skilled in the art how the present invention can 
15 be used when more than two channels are received. Similarly, it would be appreciated 
by one skilled in the art how the present invention can be used when the symbols can 
take on one of M possible values. 

In metric combining, the weighting matrix Q is a diagonal matrix. The diagonal 
entries of this matrix are the weighting coefficients for each channel. 
20 In Interference Rejection Combining (IRC), the weighting matrix Q is the 

inverse of the correlation matrix of the impairment. Specifically, the impairment vector 
i(n) may be defined as 
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i(n) 



\(n) 



(36) 



and let R m = Efifaftfo) 11 }. In IRC, the weighting matrix Q is simply Q = r~ 1 . 
IRC is described more fully in U.S. Patent No. 5,680,419 which is incorporated by 
reference herein. 

For a multi-channel MLSE with Ungerboeck metric, consider the case with two 
channels; hence, for each "n" there would be two received samples rjn) and r b (n) as 
given by equations (28)-(29). 

In this case, the Ungerboeck branch metric is given by: 



M h [n) = Re{s*(n)[-2z(n) + ?0 ^(») + q^s^n - 1)]} 



(37) 



where z(n) is obtained from both received signals and both channel taps, and (q 0 , q x ) 
are obtained from both channel taps and the inverse of the impairments correlation 
10 matrix (see US Patent 5,680,419). It is important to note that (q 0 , q x ) do not depend on 
the hypothesized symbols and are fixed. 

From equation (37), it is evident that computing M' h (n) in the multi-channel 
case requires computing 



t h (n) = qoS h (n) + qjS h (n - 1). 



(38) 



15 Computing t h (ri) would normally require two multiplications. These multiplications can 
be avoided by pre-computing the values {s# 0 , s^} for possible symbol values s and 
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storing them in a table in memory. Any hypothetical t h can then be computed by 
adding the appropriate entries from this pre-computed table. 

It would be appreciated by one skilled in the art how the present invention can 
be used when more than two channels are received. Similarly, it would be appreciated 
by one skilled in the art how the present invention can be used when the symbols can 
take on one of M possible values. 

For the case of fractionally-spaced MLSE equalization, a (772) fractionally- 
spaced MLSE equalizer will receive two samples in each sampling interval of length T. 
Let rjn) represent the even samples of the received signal, and let r b (n) represent the 
odd samples of the received signal. Note the rjn) and r b (n) are symbol-spaced. In a 
particular implementation of the fractionally-spaced MLSE equalizer, rjn) and r b (n) are 
treated as two separate received signals (resulting from the same transmitted symbol 
stream). According to this exemplary embodiment, the two-input MLSE equalizer of 
the last section is used to detect the transmitted symbol stream. The techniques 
disclosed above for reducing the complexity of the diversity of the MLSE equalizer are 
obviously applicable to this implementation of the fractionally-spaced MLSE equalizer. 

It would be appreciated by one skilled in the art how the present invention can 
be used with other forms of fractionally-spaced equalization. For example, in one 
particular realization of the fractionally-spaced equalizer, a pre-whitening filter is 
applied to each sub-sampled sequence (see Hamied and Stuber, "A Fractionally Spaced 
MLSE Receiver," IEEE 1995). This particular formulation is very similar to metric 
combining. In yet another formulation of the fractionally-spaced equalizer, an 
Ungerboeck metric is used; hence, the present invention can again be used to reduce 
the complexity of this equalizer. It would also be appreciated by one skilled in the art 
how the present invention can be used when more than two samples are received per 
symbol. 
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For a MLSE with a time-varying channel, at each iteration, a different set of 
channel taps might be considered. With the Euclidean metric, the hypothesized 
received signal is expressed as: 

r h (n) = c 0 (n) s h (n) + c x (n) s h (n - 1) (39) 

In this case, at iteration V , the table might be updated and used for computing 
all the different hypothesized received values at this iteration. A new table is computed 
for iteration "n + 1" and so forth. 

It would be appreciated by one skilled in the art how the present invention can 
be used when the Ungerboeck metric is used with a time-varying channel. The table 
would be updated every time the q/s are updated. 

For demodulation of signals modulated using offset M-ary modulation, the 
present invention can be modified to reduce the complexity of the demodulator even 
further. In offset M-ary modulation, the transmitted symbols for even "rc" are purely 
real, and the transmitted symbols for odd n n n are purely imaginary: 

s(n) = j n B t ie {1,...,M}, (40) 

where each B, is purely real. Hence, s h {n) can take one of 2M values {B^B.}^ . 

However, by pre-rotating the received signal r(n) by j" n we can reduce the 
demodulation problem to demodulating purely real symbols where each real symbol is 
in the set {B}f sl . With this pre-rotation, we would only need to store the product of 
the channels taps and each of {B }*f \s. 

It will be apparent to one of ordinary skill in the art that the present invention 
may be practiced in other embodiments that depart from these specific details. In other 
instances, detailed description of well-known methods, devices, and circuits are omitted 
so as not to obscure the description of the present invention with unnecessary details. 
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Many variants and combinations of the techniques taught above may be devised 
a person skilled in the art without departing from the spirit or scope of the invention 
described by the following claims. 
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WHAT IS CLAIMED IS: 

1 . A method for determining a branch metric in a maximum-likelihood- 
sequence-estimation equalizer which receives at least one antenna signal modulated 
with M-ary modulation, said method comprising the steps of: 

5 pre-computing values equal to a product of a complex number and a 

hypothetical symbol value; 

storing said pre-computed values in a product table; 

adding select pre-computed values from said product table to produce a 

result; and 

10 determining said branch metric using said result. 

2. The method of claim 1 wherein said complex number corresponds to a 
channel coefficient. 

3. The method of claim 1 wherein said complex number corresponds to a 
s-parameter. 

15 4. A filter in a maximum-likelihood-sequence-estimation equalizer, which 

demodulates at least one received radio signal modulated with M-ary modulation, for 
producing a hypothesized received signal sample to be used for determining a branch 
metric, said filter comprising: 

a memory for storing a product table having pre-computed values equal 

20 to a product of a channel tap estimate and a hypothetical symbol value for different 
iterations; and 

an adder for adding select entries from the product table to produce a 
hypothesized received signal sample. 
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5. A filter in a maximum-likelihood-sequence-estimation equalizer for 
M-ary modulation, said filter comprising: 

means for pre-computing a plurality of possible values to be used in 
determining a branch metric; 

a memory for storing said plurality of pre-computed possible values; and 
means for combining select pre-computed values from said memory. 

6. The filter according to claim 5 wherein said branch metric is an 
Ungerboeck branch metric. 

7. The filter according to claim 5 wherein said branch metric is an 
Euclidean branch metric. 

8. The filter according to claim 5 wherein said branch metric is a partial 
Ungerboeck branch metric. 

9. A method for computing a branch metric in maximum-likelihood- 
sequence-estimation equalizer which demodulates M-ary modulated signals, said 
method comprising the steps of: 

pre-computing a plurality of possible values to be used in the branch 
metric computation; 

storing said plurality of pre-computed possible values in a memory; 
adding select pre-computed values from said memory; and 
computing said branch metric using said added select pre-computed 

values. 

10. The method according to claim 9 wherein said branch metric is an 
Ungerboeck branch metric. 
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11. The method according to claim 9 wherein said branch metric is a partial 
Ungerboeck branch metric. 

12. The method according to claim 9 wherein said branch metric is an 
Euclidean branch metric. 

13. A method for computing a branch metric in a multi-channel maximum- 
likelihood-sequence-estimation (MLSE) equalizer which demodulates M-ary modulated 
signals, said method comprising the steps of: 

pre-computing a plurality of possible values for each channel in said 
multi-channel MLSE to be used in the branch metric computation; 

storing said plurality of possible values for each channel in separate 

product tables; 

adding select values from said separate product tables; and 
computing said branch metric using said added select values. 
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ABSTRACT 

The present invention involves demodulation of radio signals modulated with 
M-ary modulation in the presence of intersymbol interference distortion. The invention 
presents a method for reducing the number of multiplications needed to implement a 
maximum-likelihood-sequence-estimation (MLSE) equalizer for signals modulated with 
M-ary modulation. In exemplary embodiments of the present invention, the number of 
multiplications is reduced by pre-computing certain values needed for the determination 
of the branch metric and storing these pre-computed values in a product table. When a 
branch metric computation is to be made, whether it is an Euclidean branch metric 
computation or an Ungerboeck branch metric computation, certain multiplication 
operations are replaced by simple table look-up operations. As a result, the power 
efficiency and speed of the system are increased. Any receiver that demodulates 
signals that are modulated with M-ary modulation can be implemented using the 
present invention. The resulting demodulator will have a lower complexity than 
existing demodulators. 
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